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Light-cone coordinates and supersymmetric discrete light-cone quantization 
are used to analyze the thermodynamics of two-dimensional supersymmetric 
quantum chromodynamics with a Chern— Simons term in the large- A^c approx- 
imation. This requires estimation of the entire spectrum of the theory, which 
is done with a new algorithm based on Lanczos iterations. Although this work 
is still in progress, some preliminary results are presented. 
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1. Introduction 

Recent work^ has shown that thermodynamic properties can be computed 
for large- A'c supersymmetric theories. The approach is based on Ught-cone 
coordinates^ and the numerical technique of supersymmetric discrete light- 
cone quantization (SDLCQ)."^'^ Here we consider two-dimensional super- 
symmetric quantum chromodynamics with a Chern-Simons term (SQCD- 
CS),^ dimensionally reduced from three dimensions. 

Light-cone coordinates^ are defined by the time variable, a;+ = (t + 
z)/^/2, and spatial components, x = {x~,x±), where x~ = {t — z)/^/2 
and x± = {x,y). The light-cone energy and momentum are given by p~ = 
{E - Pz)/V2 and p={p+ = {E-\- Pz)/V2, p± = {p^,py)), respectively. 

For field theories quantized in terms of these coordinates, the standard 
numerical technique is discrete light-cone quantization (DLCQ).^'^ Space 
is restricted to a light-cone box —L < x^ < L, —Lj_ < x,y < L± with 
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periodic boundary conditions. Momentum is then discretized as p'^ r'^i? 
p^j_ {-^riix, -^niy) with rij, riix and riiy all integers. The limit L 
oo is exchanged for a limit in terms of the integer (harmonic) resolution^ 
K = ^P~^ for fixed total momentum P+. Because the rij are positive, the 
number of particles is limited to no more than K. Integrals are replaced by 
discrete sums. 

SDLCQ^'^ is a special form of DLCQ that preserves at least part of the 
supersymmetry algebra 

{Q+,Q+} = 2V2P+, {Q-,Q-} = 2V2P-, {Q+,Q-} = -4P^. (1) 

Instead of discretizing the Hamiltonian P~ directly, the supercharge Q~ is 
discretized, and P~ is computed from the algebra as 

-PsDLCQ = {Q 'Q } 7^ -^DLCQ- (2) 

For ordinary DLCQ, one recovers supersymmetry only in the infinite reso- 
lution limit. 

After a brief summary of the SQCD-CS theory, we show how thermo- 
dynamic quantities can be constructed from the partition function. This 
requires knowledge of the spectrum of the theory, which we obtain numer- 
ically with an iterative Lanczos algorithm. Some preliminary results are 
presented and future work discussed. 



2. Supersymmetric QCD 

We consider a dimensional reduction from 2+1 to 1-1-1 dimensions oi M = 
1 supersymmetric quantum chromo dynamics with a Chern-Simons term 
(SQCD-CS). 5 The action is 



-g [*Ae + ^U*] + -Ar'^I^^A 



(3) 



-I-kAA 



where the adjoint fields are the gauge boson A^^ (gluons) and a Majo- 
rana fermion A (gluinos) and the fundamental fields are a Dirac fermion 

^ (quarks) and a complex scalar ^ (squarks) . The Chern- Simons coupling, 
K, has the effect of providing a mass for the adjoint fields. The covariant 
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derivatives are 

D^K = d^K + ig[A^, A] , D^^ = d^^ + igA^^ , 

D^^ = d^^ + igA^^. (4) 

The fields transform according to 

5A^ = yr^A, SA=^F^,T^^■'s, 

5^ = ^e^, = -ir^sD^C- (5) 

Wc reduce to 1+1 dimensions by assuming the fields to be independent 
of the transverse coordinate x. Wc define fermion components by 

A=(a,a)^, *=(v,Vi)^, Q = {Q+,Q-f. (6) 
There are constraints, which in hght-cone gauge (A+ = 0) are written 

a_A = -^([A^A]+^eV^-^V^t), (7) 

a_Vi = -^AV+^A€-«A/v^, dlA-=gJ, (8) 

with 

J = d-A'^] + ^{A, A} + Kd-A^ - ihd-(,i^ + i^d-C^ + s/2i)ip^ . (9) 
The reduced supercliarge is 

Q~ =9 j dx- |23/4 (^[A\d-A^\ - Kd-A^ + ^{A,A}^ ^A 

--^ (iV2^d-£} - iy/2d-^£} + 2^^-^) (10) 
-2{i^A^^ + i,^A^i)] . 

In the large- iVc approximation, there are only single-trace Fock states, 
the mesons 

fl{ki)al,,{k^) . ..bl,^^^{K-i) . ..fliknm (11) 

and glueballs 

Tr[a^(A:i)...6|„,„^^(A:„)]|0), (12) 

where fj and // create fundamental partons and aj^ and fej^ create adjoint 
partons. Either type of state could be a boson or a fermion. 
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The theory possesses a useful Zi symmetry^ 

aij{k,n^) ^ -aji{k,n^), hij{k,n^) ^ -h^i{k,n^), (13) 

which further divides the Fock space between states with even and odd 
numbers of gluons. We then diagonalize in each sector separately. 



3. Finite Temperature 

From the partition function, Z = e~P°^'^, we compute the bosonic free 

energy 

and the fermionic free energy 

:F, = -^±r dpo^4^ m (l + e--/-) . (15) 

The total free energy, once expanded in logarithms and po integrals are 
performed, is given by 

The sum over I is well approximated by the first few terms. We represent 
the sum over n as an integral over a density of states p: J p{M)dM. 

The density is approximated by a continuous function, and the integral 
/ dM is computed by standard numerical techniques. 



4. Lanczos Algorithm for Density of States 

The discrete density of states is p(M^) = X^„c?„(5(M^ - M^), where d„ is 
the degeneracy of the mass eigenvalue M„. The density can be written in 
the form of a trace over the evolution operator e~ 



1 

p^M^) = —L^ / e*^'^^/2^+TVe-'^"^^rfa;+. (17) 

We approximate the trace as an average over a random sample of vectors^ 

1 ^ 

p{M^)^-Y,Ps{M\ (18) 
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with ps a local density for a single vector |s), defined by 
1 

"'^^'^^4^ J e'^''^/2^^(s|e-^^"-^|s)rfa;+. (19) 

The sample vectors \s) can be chosen as random phase vectors;^" the co- 
efficient of each Fock state in the basis is a random number of modulus 
one. 

Wc approximate the matrix element (s|e^'^ ^ \s) by Lanczos 
iterations. Let D be the length of \s), and define = -^\s) as the 

initial Lanczos vector. The matrix element (ui|e~*^ can be approx- 

imated by the (1,1) element of the exponentiation of the Lanczos tridiag- 
onaHzation of P~ . Let P~ be the tridiagonal Lanczos matrix. It can be 
exponentiated by first diagonalizing it: 

n-cj = (20) 

such that = UAU-\ with Uij = (cpi and Aij = 6ij^. The (1,1) 
element is given by 

= E l«)lpe-^^^n^"/^^^ (21) 

n 

The local density can now be estimated by 

P,(M2) ~ ^snSiM^ - Ml), (22) 

n 

where Wgn = -D|(c^)ip is the weight of each Lanczos eigenvalue. Only the 
extreme Lanczos eigenvalues are good approximations to eigenvalues of the 
original . The other Lanczos eigenvalues provide a smeared representa- 
tion of the full spectrum. 

Prom this density of states, we compute the cumulative distribution 
function (CDF), N{M^) = dM^p{M^) as an average 

iV(M2)::, i^iV,(M2), (23) 

S 

of local CDFs 

dM^Ps{M^) ~ J^WsnOiM^ - M^). (24) 



The convergence of the approximation is dependent on the number of 
Lanczos iterations per sample, as well as the number S of samples. Test 
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runs indicate that taking 20 samples is sufficient. The number of Lanczos 
iterations needs to be on the order of 1000 per sample; using only 100 leaves 
errors on the order of 1-2%. 

5. Preliminary Results 

Some preliminary results are presented in the accompanying figures. Figure 
1 compares the numerical results for the CDF to analytic results, which can 
be obtained when the Yang-Mills coupling is zero. The numerical results 
are quite good, with only one noticeable deviation, at large M^, a region 
where the discrete spectrum is sparse. Figure 2 shows the free energy at 
particular values of the Yang-Mills coupling. 
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Fig. 1. Cumulative distribution function for resolution ii" = 13 in the analytically 
solvable case of zero Yang-Mills coupling. The numerical and analytic solutions are 
compared. 



6. Future Work 

Additional work is in progress to complete this study of finite temperature 
properties of two-dimensional SQCD. Beyond this particular effort, one can 
consider finite- A^c effects, with baryons and mixing of mesons and glueballs, 
and the full three-dimensional theory. As these techniques mature, analysis 
of four-dimensional theories can be considered. 
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Fig. 2. Free energy at fixed Yang-Mills coupling g as a function of temperature T for 
resolution K = 14. 
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